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Abstract 

The mechanism for phase distortion in linear beam vacuum devices 
is identified in the simplest system that pertains to such devices, the 
force-free drifting electron beam. We show that the dominant cause of 
phase distortion in a force-free drifting beam is the inverse dependence 
of arrival time on velocity for an electron, i.e., the nonlinearity," 
and that a secondary cause of phase distortion is the nonlinearity of the 
velocity modulation. We claim that this is the mechanism for phase dis- 
tortion in all linear beam vacuum devices, and provide evidence from a 
traveling wave tube calculation. Finally, we discuss the force-free drifting 
beam example in the context of the "self-intermodulation" description of 
phase distortion recently described in Refs. [J. Wohlbier and J. Booske, 
Phys. Rev. E, vol. 69, 2004, 066502], [J. Wohlbier and J. Booske, IEEE 
Trans. Elec. Devices, To appear.] 

We recently reported on mechanisms for phase distortion in linear beam vac- 
uum devices [3,4]. The mechanism for phase distortion was reported there as 
a "self-intermodulation" process, where harmonic beam distortions interacted 
with the fundamental beam modulation to produce a phase shift at the funda- 
mental. While this frequency domain view of phase distortion is a useful one, 
we felt that a corresponding time domain view of phase distortion would be a 
useful contribution to the overall understanding of phase distortion. In this note 
we consider the simplest possible system that pertains to linear beam electron 
devices, a space charge free ( "force- free" ) drifting electron beam. Indeed the 
physical mechanism for phase distortion exists in the force-free drifting beam, 
and hence the system provides the most lucid example in which to study phase 
distortion. 

A 1-d force-free drifting beam is described by Burger's equation for the beam 
velocity u, 

u t + uu x = 0, (1) 

where x is space, t is time, and the subscripts indicate partial derivatives. To 
apply Eq. (JTJ to a klystron where the space charge force is negligible, for exam- 
ple, one sets the boundary value of u as the sum of a dc beam velocity uq and 
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a sinusoidal perturbation due to the cavity modulation. Assuming that the dc 
velocity is normalized to 1, this is written 



u(Q,t) = 1 + esinwi. (2) 

Given the solution to Eq. with boundary condition @, the beam density 
evolution is obtained by solving the continuity equation 

Pt + {pu) x = (3) 

for an appropriate density boundary condition. 

The force-free Burger's equation implies that the velocity of an electron 
(fluid element) does not change, i.e., the (x,t) trajectories, or characteristics, 
of Eq. Q are straight lines with slopes determined by the boundary data. In 
vacuum electronics the density bunching that results from the sinusoidal velocity 
modulation is called "ballistic bunching." Initially we will consider inputs for 
which there is no electron overtaking, and inputs such that electron overtaking 
occurs will be considered later. 

In principle the force-free drifting electron beam can be solved exactly prior 
to electron overtaking, although closed form analytic solutions need to be writ- 
ten in terms of infinite series [1]. Equation is solved using the method of 
characteristics. The method of characteristics involves changing the equations to 
Lagrangian independent coordinates (material coordinates), solving the equa- 
tions in Lagrangian coordinates, and changing the solution back to Eulerian 
coordinates. For a fluid element that crosses x = at time to we write the 
transformation from Lagrangian to Eulerian coordinates as the time function 
t(x,to), i.e., the time fluid element to arrives at x, with i(0,io) = to. Since the 
velocity is a constant for each fluid element, the velocity solution in Lagrangian 
coordinates is 

u{x,to) = 1 + esinwto. (4) 

That is, independent of location x, the velocity of an electron is set by the 
time to at which it crosses x — 0. Since the electron orbits are straight line 
trajectories in (x,t) space, we can infer the solution for the function t(x,to) 

t(x,t ) = h+ . X + . (5) 
u(x,t ) 

= to + — — -■ (6) 

1 + e sin Loto 

To solve for u(x, t) one needs to invert the function t(x, to), i.e., calculate to(x, t), 
and substitute it into Eq. Since we restrict, for now, our attention to 

input levels such that no electron overtaking occurs, i.e., the characteristics 
do not cross, Eq. © can be inverted. Unfortunately t(x, to) in Eq. is a 
transcendental function of to, so an analytic inverse needs to be expressed in 
terms of an infinite series. 
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For phase distortion we are ultimately interested in the beam current and 
beam density, since they linearly drive circuit or cavity fields. The continuity 
equation in Lagrangian coordinates is given by [1,5] 



p(x,t ) 



p(0,t o )u(0,t o ) ^ 



dt 



u(x, t a ) 

p(0,*o) 



at 

dt 



(8) 



where the last equality comes from the fact that w(0,io) = u(x,to) since we 
are considering a force-free drifting beam. For the density solution in Eulerian 
coordinates one composes the density in Lagrangian coordinates JHJ) with the 
mapping from Eulerian to Lagrangian coordinates to(x, t). 

The factor dt/dto is the Jacobian of the transformation from Lagrangian to 
Eulerian coordinates, and it quantifies the amount of stretching or compressing 
in time the electron beam undergoes. That is, for a fluid element entering the 
system of time length Sto, it will occupy a time length of St = (dt(x, to) / dto)5to 
downstream at position x. For purposes of explaining phase distortion it will 
be convenient to write the Jacobian in terms of derivatives with respect to both 
to and u, i.e., 1 

dt dt dt du 

dt dto du dt ' 

In Fig. ^ we show density solutions versus time, where a root finder was used 
to invert Eq. The results are for the velocity modulation @ with uj = 2n 
and a uniform density boundary condition p(0, to) — 1. To accentuate the effect 
of phase distortion we use a strong velocity modulation, e = 0.389, and observe 
the output at x = 0.3. There is no electron overtaking at x if the Jacobian stays 
greater than zero, i.e., when 

ewcoswto . * , in s 

X ; j ^ 1 (10) 

(1 + e sinw^o) 

for all to. In the present case the maximum value of the left hand side of Eq. l|lUf> 
is equal to 0.9992 with x — 0.3. The linear velocity and density solutions are 
solutions to the (normalized) linear equations 

u t +u x = (11) 
Pt + Px+Ux = (12) 



1 We use the function i(x,to,u), cf. Eq. JSJ, to formally differentiate from the function 
t(x,to), cf. Eq. (|>J, so that the left hand side of Eq. J^J and the first term on the right hand 
side of Eq. are not considered equal. We did not use the hat ( " ) notation in Eqs. and 
151 where there was no chance of confusion, and we will drop the notation for the remainder 
of the paper. 
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Figure 1: Nonlinear, linear, and fundamental component of nonlinear density 
solution versus t at x — 0.3 for two periods. 

and, for the same boundary data as used for the nonlinear problem, are given 

by 

U\i n (t,x) = 1 + e sinoj(t — x) (13) 
pn n (t,x) = 1 + xeu cos u>(t — x). (14) 

In Fig. ^ the nonlinear density solution, the linear density solution, and the 
fundamental component of the nonlinear density solution obtained by using a 
fast Fourier transform are shown. The phase shift between the linear density 
solution and the fundamental component of the nonlinear density solution is the 
"phase distortion" of the nonlinear solution. It is clear from the figure that the 
reason for the phase shift of the fundamental component of the nonlinear density 
solution from the linear solution is because the density is larger "on average" 
behind (in time) the density peak (early times are to the left). Therefore, to 
explain phase distortion one must explain the reason for the asymmetry in the 
density about the density peak. (Also note that for this value of e the nonlinear 
solution shown in Fig. ^ shows "amplitude distortion" in that the fundamental 
component of the density is smaller than the predicted linear density amplitude.) 

The reason for the the higher density behind the peak comes predominantly 
from the "1/u nonlinearity" in the arrival time function and secondarily 
from the nonlinearity of the esinwio velocity modulation. It is intuitive that 
about a larger initial velocity u, a given deviation ±Su will result in a smaller 
=p<5i at x than the same ±Su would have on =p<5i for a smaller initial velocity u 
[a fluid element with velocity u ± 5u arrives earlier (later), at time t =F St, than 
a fluid clement with velocity u which arrives at t]. This intuition is of course 
borne out by the derivative 



evaluated at large and small values of velocity u. This dependence causes phase 
distortion in that relatively slower fluid elements stretch and compress differently 
than the relatively faster fluid elements. That is, the size of density fluctuations 
at a fluid element depends on its initial velocity, with smaller initial velocities 
having larger fluctuations. This effect is best demonstrated with an example 
where more specific points regarding the interplay of the 1/u nonlinearity and 
the nonlinearity of the velocity modulation may be emphasized. 

In Fig. [3 we show the velocity modulation, the Jacobian, and the density 
versus fluid element t at x — 0.3. Note that since u(x,t ) — u(0,t o ) the 
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Figure 2: Beam velocity, beam density, and Jacobian versus fluid element la- 
bel to at x = 0.3. That is, the value of velocity, density and Jacobian that the 
fluid element that crossed x = at time to will have when it arrives at x = 0.3. 
Results show slightly more than one period in to. 

velocity versus fluid element label to is also the initial velocity modulation at 
x = versus fluid element label to- In the region B-D the beam is bunched 
because fluid elements entering at any time are given a faster initial velocity 
than the fluid elements entering just prior to them, where again we assume 
for the time being that modulations are not strong enough to cause electron 
overtaking. It is instructive to look at the expression for the Jacobian, 

dt dt du . . 

oto du ato 

and determine, for example, the fluid element at which the density will be 
maximum. The Jacobian is inversely proportional to density, and hence the 
maximum density occurs the minimum Jacobian; for a given Stg the minimum 
5t corresponds to the maximum fluid compression. Consider the two derivatives 
on the right hand side of Eq. I|16fl. As described above, due to the inverse 
dependence of t on it, for a given 5u fluid elements with relatively larger initial 
velocities in B-D will be compressed less than fluid elements with relatively 
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smaller initial velocities. That is, dt/du is largest negative at B and smallest 
negative at D. Furthermore, for the sinusoidal velocity modulation the change in 
velocity for a given <5to, du/dto, is zero at point B and increases to a maximum 
at point C. The minimum Jacobian will occur at a fluid element between points 
B and C such that (dt/du)(du/dto) is the largest negative. The value of t for 
which this happens can be computed by setting the derivative of the Jacobian 
with respect to to to zero, and solving for to. For the linear solution point C is 
the fluid element around which the density is maximum. 

The density asymmetry about the peak also comes from the 1 ju nonlinearity. 
Since the compression is enhanced in B-C relative to C-D, region B-C becomes 
smaller than C-D, as seen in Fig. EI The results in Fig. [31 are the solutions in 
Eulerian coordinates which are obtained by composing the Lagrangian solutions 
in Fig. [2] with the map to(x,t), i.e., the mapping that determines which fluid 
element to arrives at a given (x, t), shown in Fig. 01 Indeed while the maximum 
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Figure 3: Beam velocity, density, and Jacobian versus t at x = 0.3 for one period 
of t. Results are those of Fig. [21 composed with to(x, t). shown in Fig. 0] 

density fluid element lies nearly half way between points B and C at the input 
(see Fig. |2f ) it is very near to point C by x = 0.3, the point about which 
the linear density is maximum, as seen in Fig. [21 Because B-C is much more 
compressed than C-D, loosely speaking, the regions that neighbor the region of 
maximum density are a stretched region A-B earlier in time and the compressed 
region C-D along with the stretched region D-A (periodic waveform) later in 
time. Furthermore, the stretched region A-B has enhanced stretching over 
region D-A, contributing to the higher density behind (in time) the peak. 

In sum, the above describes how phase distortion of a force- free drifting beam 
is a manifestation primarily of the inverse dependence of the arrival time t(x,to) 

2 Since the velocity of a fluid element does not change in the force- free drifting beam 
u(x, to) = u(0,to). Hence u(x,to) in Fig. 121 can be used to determine the phase position of 
fluid element to with respect to the velocity modulation at x = 0, even though the figure 
caption indicates that the quantities are evaluated at x = 0.3. 
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Figure 4: Mapping to(x,t) from Eulerian to Lagrangian coordinates at x = 0.3 
for one period of t. 



on velocity u, but also depends on the nonlinearity of the velocity modulation. 
To prove that the "dominant" cause for phase distortion is the 1 /u nonlinearity 
we do the following calculation. First, for a small velocity modulation (e small) 
we can linearize the 1/u nonlinearity to get 

« t + x(l - esinwt ), (17) 
w 1 — xeuj cos toto- (18) 

For this arrival time function the location of maximum density (minimum Jaco- 
bian) is the same as for the linear solution, point C in Fig.0 However, with this 
approximation the problem does not immediately reduce to the linear solution 
since the characteristics do not all have the same slope. 3 In this limit there is no 
phase distortion since Eq. (|18fl is an even function about the Jacobian minimum, 
and hence the density is an even function about the density peak. Furthermore, 
we could replace the modulation by any periodic function that is odd in to (plus 
or minus an arbitrary time shift) for a period and get the same result. Thus, 
by linearizing the 1/u nonlinearity and not linearizing the velocity modulation 
the phase distortion is removed. This confirms that the 1/u nonlinearity is the 
dominant cause of phase distortion. 4 

Up to now we have restricted the input to a level such that electron trajecto- 
ries do not cross. The reason for such a restriction was to ensure that the phase 
distortion physics was not potentially clouded by electron overtaking physics. 

3 To get to the solution of the linear problem one can substitute Eq. 1181 into Eq. JSJ, use 
that e is small to move the denominator to the numerator with a sign change on the cosine 
term, and use to = t — x which can be obtained by taking e = in Eq. 1171 . 

4 One might be tempted to linearize the modulation without linearizing the 1/u nonlinearity 
for a complementary view. However, it is not possible to linearize the modulation and maintain 
its periodicity since a periodic function is nonlinear in its argument. 



t{X,to) 

dt_ 
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In fact, no such restriction was necessary, and the principle of phase distortion 
is the same even when the input is set such that electron trajectories cross. In 
Fig. Owe show the beam velocity, nonlinear beam density, linear beam density, 




Figure 5: Nonlinear beam density, linear beam density, fundamental component 
of the nonlinear beam density, and nonlinear beam velocity for an input such 
that electron overtaking occurs. 

and the fundamental component of the nonlinear beam density for an input 
of e = 0.6. The velocity solution and the two peaked structure of the density 
confirm the multi- valued nature of the solutions. From the fluid element labels 
we see that region B-D is in the multi-phase region of the density, and that 
the higher density behind the peak is again because region D-A (or region B-A 
outside the multi-phase zone) is not stretched to the extent that region A-B (or 
region A-C outside the multi-phase zone) is. 

Even though the preceding analysis considered only a force-free drifting 
beam, we claim that phase distortion mechanism identified holds more gen- 
erally for klystrons when space charge forces are considered, and in traveling 
wave tubes (TWTs) . The reason is again the inverse dependence of arrival time 
on velocity. That is, electrons slowed down from the dc beam velocity will have 
enhanced stretching or compression over those that are sped up from the dc 
beam velocity. Although we do not provide an exhaustive analysis, for illus- 
tration we consider the electron beam density versus time from a Lagrangian 
TWT model [6] in Fig. EJ Even though electrons experience forces due to space 
charge and circuit fields in a TWT, as can be inferred by comparing the beam 
velocities of Figs. Eland El the cause of phase distortion is the same. Electrons 
slowed down from the dc beam velocity will spread more in time than electrons 
sped up from the dc beam velocity. The result is that the density wave form 
has a relatively higher density behind the peak than in front of it, as seen in 

Fig. ED 

In Rcfs. [3,4] we described phase distortion from a frequency domain per- 
spective as a "self-intermodulation" process, whereby harmonic distortions mix 
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Figure 6: Beam density and beam velocity versus time for two periods from 
a Lagrangian TWT calculation accounting for space charge and circuit fields. 
The expanded view of the beam density is to highlight the density asymmetry 
about the peak. 



with the fundamental to produce distortions at the fundamental. Below we 
outline this view of phase distortion so that it may be compared to the time 
domain view given above. 

If we express the velocity u and density p with Fourier series (for periodic 
inputs), then the fundamental frequency components of nonlinear products of 
p and u are seen to come from mixing of the second harmonic with the funda- 
mental frequency, and mixing of higher order harmonics. In particular, if we 
define state variable envelopes as in Refs. [3,4], e.g., 

oo 

u(x,t) = J2 fimWe*™ 1 *-*, (19) 

rn — — oo 

then the continuity equation for the force-free drifting beam gives 

iluiu n + i muju m u n 

- Uuj ^2 u m p n (20) 

m + n=t 

where the approximations used in Refs. [3,4] have been made. From Eq. I|20[l 
products of frequencies such as (m, n) = (—1, 2), (2, —1), (—2, 3), (3, —2), etc. are 
seen to influence the fundamental frequency In Refs. [3,4] we considered 

(m,n) = (— 1, 2), (2, — 1) (third order intermodulation, "3IM") and (m,n) — 
(—2, 3), (3, —2) (5IM) contributions to phase distortion in linear beam devices. 

Equation (|20|l together with the corresponding nonlinear envelope equation 
for the velocity has an analytic solution that may be expressed in an infinite 
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series of complex exponentials. The first terms in the series correspond to the 
linear solution as seen in Fig. The linear terms are used in the equation for 
the second harmonic envelopes to produce the first terms in the series solution 
at the second harmonic. The second harmonic terms are then combined with 
the linear terms from the fundamental solution to produce the next set of terms 
at the fundamental. These complex exponentials add to the linear density solu- 
tion and produce a phase shift (distortion) in the density solution. As it turns 
out, this process of generating harmonics and then additional terms at the fun- 
damental must be continued to get an accurate representation of the nonlinear 
phase shifted density seen in Fig. ^ F° r the example in Fig. ^ we found that 
terms higher than 11 IM were required to adequately approximate the nonlinear 
density, whereas in Refs. [3,4] we found that 5IMs were sufficient to predict the 
phase distortion. The difference between the two cases is that the case in Fig.^ 
has a much larger relative input modulation. 

We have identified the mechanism for phase distortion in the simplest sys- 
tem that pertains to linear beam vacuum devices, the force-free drifting electron 
beam. The dominant cause of phase distortion is shown to be the inverse depen- 
dence of arrival time for an electron on its velocity, i.e., the "1/u nonlinearity," 
and that a secondary cause of phase distortion is the nonlinearity of the veloc- 
ity modulation. Although we show this to be the case for the force-free drifting 
beam, we claim that the 1/u nonlinearity is also the cause of phase distortion in 
other linear beam devices, such as the traveling wave tube. That is, even though 
electrons are experiencing forces from the circuit and space charge fields, the fact 
remains that charges slowed down from the dc beam velocity will stretch and 
compress more than charges sped up from the dc beam velocity. Results from 
a traveling wave tube calculation are given that show the same characteristic 
asymmetric density modulation that is seen in the force-free drifting beam. We 
also show that the 1/u nonlinearity is the cause of phase distortion regardless 
of whether drive levels are strong enough such that electron overtaking occurs. 

The identification of the mechanism for phase distortion suggests how inputs 
might be tailored to ameliorate phase distortion. The obvious candidates would 
involve somehow reducing the amplitude of the velocity modulation on the neg- 
ative half cycle to lessen the effect of the "enhanced stretching," or to provide 
a density modulation 180° out of phase from the velocity modulation so that 
the "enhanced stretching" before the peak starts from a higher density value 
relative to the density behind the peak, and would be balanced about the den- 
sity peak at the output. The latter scheme may be facilitated by cold cathode 
technology [2]. It is anticipated that either of these schemes, and potentially 
any scheme, may come with an associated gain compression, as in Ref. [3] where 
harmonic injection was used to set the fundamental output phase at a cost of 
reducing output power. It is also possible that further study may show that any 
tailoring of the velocity modulation may just be a manifestation of harmonic 
injection, proving the usefulness of the complementary views of phase distortion 
given in this note and in Refs. [3,4]. 
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